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Abstract 

Vibrations in modem structural and mechanical 
systems can be reduced in amplitude by increasing 
stiffness, redistributing stiffness and mass, and/or adding 
damping if design techniques arc available to do so. Linear 
Quadratic Regulator (or LQR) theory in modern 
multivariable control design, attacks the general dissipative 
elastic system design problem in a global formulation. 
The optimal design, however, allows electronic 
connections and phase relations which are not physically 
practical or possible in passive structural-mechanical 
devices. This paper addresses the restriction of LQR 
solutions (to the Algebraic Riccati Equation) to design 
spaces which can be implemented as passive structural 
members and/or dampers. A general closed-form solution 
to the optimal free-decay control problem is presented 
which is tailored for structural-mechanical systems. The 
solution includes, as subsets, special cases such as the 
Rayleigh Dissipation Function and total energy. 
Weighting matrix selection is a constrained choice among 
several parameters to obtain desired physical relationships. 
The closed-form solution is also applicable to active 
control design for systems where perfect, collocated 
actuator-sensor pairs exist . Some examples of simple 
spring mass systems are shown to illustrate key points. 

L Introduction 

Many modem structural and mechanical systems 
are subject to dynamic input forces which excite flexible 
vibrations. These vibations can be reduced in amplitude 
by increasing stiffness. However, reducing mass or inertia 
at the risk of reducing stiffness is an important factor for 
other reasons, such as launch costs or rigid-body control 
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energy. A third design parameter, damping, offers an 
alternative for the stiffness-versus-mass dilemma, if 
suitable design methods can be developed for dissipative 
structures. Also, redistribution of stiffness and mass 
without increasing total mass can be beneficial if design 
techniques are available to do so. 

Linear Quadratic Regulator (or LQR) theory in 
modem multivariable control design, attacks the general 
dissipative system design problem in a global formulation 
based first-order equations. This formulation usually 
yields a requirement for computer solution of a nonlinear 
matrix equation (the Algebraic Riccati Equation). For the 
case of structural and mechanical systems where second- 
order governing equations predominate, the control 
equations are twice the usual matrix dimension, dependent 
on arbitrarily selected weighting matrices, and are difficult 
to interpret in terms of physical design parameters. In 
control design, a computer can implement connections and 
phase relations which are not physically practical or 
possible in the structural-mechanical case (e.g., Hanks and 
Skelton 1 , and Hvorat 2 ). Nevertheless, the theory is 
applicable to passive structural systems and the challenge 
in using it is to restrict it to physically possible design 
spaces. Previous studies, such as Skelton, et. al. 3 *^, 
have considered physically constrained computer solutions. 
Others, such as Kojitani and Ikeda 6 and Belvin and Park 7 
have provided closed-form solutions for special cases. 

This paper presents a general closed-form solution 
to the optimal free-decay control problem which is 
constrained to predetermined structural-mechanical 
geometry. It is derived entirely in physical coordinates 
using mass, stiffness, and viscous damping matrices. The 
solution results in a final (closed-loop) objective function 
which is a generalized energy expression for dissipative 
structural-mechanical systems. Application is discussed to 
special cases including the Rayleigh Dissipation Function, 
total energy, and known closed-form solutions for 
particular objective functions. Weighting matrix selection 
is a choice among several parameter matrices. For design 
objectives commonly used in LQR control system design 
(i.e., choices of Q weighting matrices), direct solution for 
optimal design parameters (R weighting matrices) is 
possible, but frequently structurally impractical. The 
closed-form solution and its physical interpretation is also 
applicable to the active control problem for similar 
systems where perfect, collocated actuator-sensor pairs 
exist and, hence, is applicable for initial control design and 



II. Connectivity Modeling Qf Discrete 
Physical Systems 

Any discrete physical system, or discrete 
representation of a continuous physical system, may be 
modeled in connectivity form, that is, as a three-matrix 
product SGH where each column of and each row of H 
carries connection geometry information and G is_a 
diagonal "gain” matrix . G has dimension mxm where m 
is the number of connections between nodes, *B has 
dimension nxm where n is the number of input (forced) 
nodes and H has dimension mxk where k is the number of 
output (sensed) nodes. In passive structural/mechanical 
systems, or active and mixed active-passive systems with 
perfect collocated sensors and actuators, where * 

denotes the complex conjugate transpose of $ . Thus, the 
usual nxn symmetric mass, stiffness, and damping 

parameter matrices, M, and X become 

M = <BG X. , and X = ( 1 ) 

where 

Gj = diagf-g,-], (2) 

are mxm diagonal matrices. The number of terms, g^ on 
the diagonal of G in each parameter is large, i.e., on the 
order of n+(p/2), where p is the number of non-zero off- 
diagonal terms in each matrix. It is convenient to keep $ 
the same for all three matrices and enter gj=0 for non- 
existent connections, provided there is no requirement for 

Fig. 1 Simple Floating Three-Mass Elastic System 

the inverse of G to exist. The dimension of G is usually 
determined by the number of stiffness connections in % ■ 
The columns, bj of <B contain one or two non-zero entries 
of direction cosines to ground or between nodes. The 
diagonal terms of each G matrix contains local parameter 
magnitudes between nodes such as lumped mass, spring or 
damper coefficients in spring-mass systems or, similarly, 
truss elongation stiffness in trusses. For illustration, 
consider the simple system of Fig. 1. For this system, 
the parameter matrices are as follows: 
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and G# has the same form as G^. If compactness is 
important, or if inverses of G are necessary, these matrices 
may be simplified by eliminating the zero-gain columns 
and rows and the corresponding columns of but the 
form here is useful for explanation of relative positions of 
terms. Furthermore, if the springs and dampers are active 
devices, or are in parallel with active devices, the net gain 
terms can be negative or zero. 


In this paper, we will conside r two types of 
gains, the initial system gainsJ T e. M f X. and ff X .and 
"passive control” gains ( changes in mass, A^,stiffnes s, 
AX> or damping, A'Dy A control gain may involve 
changing an existing gain value or introducing a new 
connection which was not present in the initial design. 
The term closed-loop gains, denoted by a hatted symbol, 
will refer to the final mass, stiffness, and damping 

(i.e.,*/ = + AM, X = X + A9C, and T> = £>+ M>), 

after either passive or active changes are made. Active 
controls using acceleration, position, and rate feedback 
effectively produce changes in mass, stiffness, and 
damping, respectively (assuming negligible actuator and 
sensor dynamics). 
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Fig. 2 Three-Mass Floating Elastic System with Possible 
Collocated Control Connections 



In designing a set of springs (passive position- 
feedback) and dampers (passive rate-feedback), a set of 
active position, rate and acceleration feedback controls, or a 
combination of all, a large number of connections are 
possible. For the example of Fig. 1, the four possible 
additional connection paths for collocated pairs are shown 
schematically in Fig. 2. There are 24 possible connection 
combinations for the system of Fig. 2, two of which are 
shown with their corresponding ‘B matrices in Fig. 3. 
Configuration 1 includes all possible (collocated) 
connections and represents the type of solution which may 
be obtained from standard, unconstrained LQR controls 
designs. In such designs, interconnection capability is 
assumed available at all nodes but, especially in the case of 
passive implementation, is rarely practical. 






formulation, but certain parameters such as weighting 
matrix values, can be physically interpreted using the 
connectivity concept. 


III. Theo ry Preliminaries 

Linear Quadratic Regulator (LQR) control theory 
has potential for structural dynamics design in mechanics. 
This theory is widely available in textbooks, such as 
Skelton^, but an outline is presented here for purposes of 
reference and terminology in deriving later results. Free- 
decay and full-state feed-back motion will be discussed. 

The time-invariant LQR approach minimizes an 
objective function, V, which is the infinite-time integral 
of the sum of weighted quadratic potential function of the 
state vector, z, and the control vector, u. It is well known 
that the minimization of 


V = — J"(z * Qz + u * Ru)dt 
2 


(3) 


subject to the constraints of the system equations of motio 
z = Az + Bu (4) 

yields 
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Configuration 1 


Configuration 2 


KA + A * K - KBR ‘BK + Q = 0 


( 6 ) 


Fig.3 Some Alternative Connection Geometries for a 
3-Mass Simple System 

Hanks and Skelton 1 , using a finite-element, 
cantilever beam example, showed that the addition of 
connections from analysis nodes to ground and between 
non-adjacent analysis nodes are required by an LQR design 
solution. Such connection freedom may not be cost- 
effective, even if it is possible. Hvorat^, studying 
minimization of squared displacements on a two-mass, 
simple system, concluded from mathematical criteria that 
an LQR controller could not be implemented by a single, 
passive, non-grounded controller. The option to constrain 
the solution to allow only physically practical connection 
locations, such as in Configurations 1 and 2 in Fig. 3, is 
necessary in practical passive or active control design. 
Skelton, et. al.^. 4 , using the above diagonal-gain 
connectivity approach to fix control configuration 
geometry, applied an iterative constrained output feedback 
solution to minimize a quadratic objective function while 
satisfying predetermined limits on parameter changes. The 
following sections discuss a closed-form solution to the 
LQR problem which provides the optimal solution while 
maintaining predefined passive, or collocated active, 
control geometry. It is not restricted to the diagonal gain 


Substituting the linear optimal feedback control (5) into 
(4), the closed-loop system equations of motion become 

z = Az-BR ‘BKz = (A + AA)z = Az (j) 

Similarly, (3) yields the optimal objective function 

V = — J„"(z * Qz + z * KBR 'BKz)dt 
2 

= -Jj"z*(Q + AQ)zdt 
2 

= — J„" (z * Q z)dt = — z * (O)Kz(O), 

2 2 

Q = Q + KBR 'B * K (g) 

where z(0) is the initial state vector. 

If (A,B) is stabilizable (that is, the uncontrollable 
modes are stable) and (A,Q) is detectable (that is, the 
unobservable modes are stable), then the feedback gain 
matrix, K, that solves (6) is positive semi-definite and 

A 

A = A + AA is also stable, Equation (6) can also be 
written 



4 


K A+ A* K + Q = 0 (9) 

From Liapunov stability theory, equation (9) represents a 

A 

stable system if K is positive definite and Q is positive 

A 

semi-definite, that is, K >0, Q > 0. 

Now, consider the familiar equation 


*• • 

*fq+(© + £)q+:*Sl-'BU (jo) 

where and Q are mass, stiffness, damping and 

gyroscopic-coupling matrices, respectively, having the 
following properties: 


M = M , >0, AC = 3C* > 0, 
© = ©* 2: 0, and Q = -Q * 


(ID 


z * Qz = z*Qz + u*Ru = z*(Q + KBR 4 BK)z 
where the 2n x 2n matrix Q reduces to 


(17) 


* fQ„ Q»1 IXx*®'' 

Q = + R (®K K *»f K 

Lq d q«J U,x®'J 


(18) 


Now, for simplification of notation, we define 


r = «R * (19) 

The 2n x 2n Riccati equation, (6), from which the control 
gains arc derived, can be written in second-order form as 
four n x n matrix equations of which three are 
independent^. These three, including the notation (19), are 
as follows: 


By writing the state z, in (4) as 


1^21 Kjj + Qjj — 0 (20) 


z = 



-K 22 M'\d + g)-(D + g)*M'K 22 
(17) -i-K j 2 + K 2 , — = 0 (21) 


where q is displacement and q is velocity, the equations of 
motion (10) become, in the first-order form of (4), 


K„ - K |2 W"(© + Q) - 3C*f "K k - K n M 'r!M 'K n + Q 12 = 0 

(22) 


z = 



(A + AA)z = Az 


IY. Closed-Form Solution s Of Algebraic 

Riccati Equations 
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(13) 


The optimal control (5) is 

r k 

u„=-R’[o 'S*^ , ] K ‘ 



= -R (K 21 q + K b q) 


(14) 


and the objective function (3) has the form 


Equation (16) provides the basic relationships 
between the physical system parameters and the 
mathematical design parameter matrices. These are as 
follows: 

AflC = T&f Kjp A{2> + £) = Ti&f K 22 (23) 

Hence, solution of (20) through (22) to obtain K 12 and 
K 22 provides the optimal system stiffness and damping 
changes for any permissible set of Q and R matrices. 
Usually, these equations are solved numerically by 
computer. However, closed-form solutions are known for 
some special cases as discussed in the following section. 


V = 





4 * u * Ru 


dt 


(15) 


Substituting the optimal control (14) into (10) gives the 
equations of motion of the closed- loop, or revised system 

q+ (© + Q + m ' *f 'K a )q 

+(3C + 'BR l 'B*^r , K 21 )q-0 (j 6 ) 

A 

Similarly, Ihe closed-loop integrand z *Qz of the optimal 
objective function (15) becomes 


The Solution of Kfliilani and Ikeda 

Kojitani and Ikeda^, propose the following 
weighting matrix: 


Q = a 


~ 2 ac + »r g 

b 

g * — (2D + bT) - I'M. 


La J (24) 

Where a and b are arbitrary scalars which must be restricted 
for existence of the Riccati equation solution. This results 
in an exact solution to the Riccati equation of the form 
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K = b 


a a 

3C + a r + -2> —M 
b b 


a 

b 


b > 0 




(25) 

(26) 

u = -R''2 * (aq + bq) = -R^2 * [ al : bl ]j • i 

UJ (27) 

Proof that (25) is a solution is a simple matter of 
substituting the terms into (20) through (22). 


where the submatrices 

K 12 = K 2 , = aflf, 

lead to the very simple control vector 


Assuming that T is preassigned, the above 
solution to (20) through (22) has two arbitrary parameters 
a and b, resulting in one added constraint, K 21 =(a/b) K22- 
This constraint requires that the stiffness change be 
proportional to the damping change. Replacing this 
constraint by the condition a =0 is also possible. ( b=0 is 
not permissible if K is to be positive definite.) It also 
provides a direct implementation of control feedback loops 
in which measurements are preweighted by a scalar 

constant rather that the matrices 94.' and/or Kjj. This 
allows observation of a very important point. The closed- 
loop system equations of motion are 


Jfq+ (© + Q + bl~)q+ (3C + aOq 

- !»rq+(® + $ + AO)q+(JC + A3C)q -0 (28) 

From (28), it is clear that the closed-loop stiffness and 
damping matrix changes are entirely a function of free 
parameters (a, b, T.). The optimal control design has 
been obtained and its form is predetermined by the scalar 
constants a and b and the matrix r (composed of the 
assumed ‘B and R). Even with the luxury of full-state 
feedback, we have no indication of how to choose the key 
parameters which define the optimal controller design (i.e., 
stiffness and damping). In fact, the following proposition 
can be stated by inspection of (26) through (28): 

Proposition: Any solution of the full-state feedback 

control problem which results in Riccati solutions of the 
form (26) transforms the problem into the symmetric form 
of the standard measurement feedback problem, with R'^ 
being the equivalent of the optimal gain sought. 

Problem : Applying the diagonal gain formulation of 
Section l, R~ l becomes a diagonal gain matrix whose 

diagonal elements if 1 are the gains of single physical 
connections. For any set of mechanical! structural 
connections, find an optimal R matrix which provides the 
gains for the optimal control problem. 


The Solution of Belvin and Park, 


In the above solution, Q contains terms of the 
optimal control which is being sought. A solution which 
does not have this characteristic was offered by Belvin and 
Park? for the special case where B is square and of full 
rank. Their closed-form solution to (20) through (22) 
assumes Q and R as follows: 



R = *8 * ?C 1 

(29) 


The solution is presented in the following form: 


1\94 


rj94 

qMl c"'F. 


(30) 


Where a and P are arbitrary scalar constants and 

jj = Vl + a-l, £ = -\/2 ti + /3, 

.1 .ni 111 

94 1 \ 94 2 = 94 i, D a 94 1 

) (31) 

((D a denotes the critical damping matrix**.) This solution 
yields a control law 

* 

®u = -TJ«Jl-^q (32) 

which minimizes the sum of the strain and kinetic energies 
of the open-loop system. It is interesting to note that for 
a single-dcgree-of-freedom system the velocity feedback 

gain is , that is, proportional to critical damping. 


and 'P = 




A General Energy-Rela ted Solution, 

Equations (20) to (22), the Riccati sub-equations, 
provide three relationships among the design parameter 
matrices. There arc seven undefined parameter matrices in 
these equations, 3Qij's, 3 Kjj’s and r. Picking any four 
defines a unique relationship among the remaining 
matrices. It is commonly assumed, as in the Kojitani- 
Ikeda and Belvin-Park solutions, that T and Qjj are 
preassigned and the equations solved for K. ( r is actually 
composed of two matrices, ‘B and R, but they never appear 
independently in the above equations and may be 
considered as a unit in the solution, with the condition that 
R'l must exist.) This is not required, however. Choosing 
physically meaningful and mathematically permissible 
forms for other combinations of four terms is possible. A 
great many logical combinations of parameters exist, some 
of which may take the form of constraint equations. 


Determination of this R matrix, will be explored in The physically simple control vector form (27) of 

Section V. the Kojitani-Ikeda solution suggests a meaningful choice 


for two of the parameters, K12 and K22. Therefore, 
substituting (27) into (20) through (22), we obtain a 
general energy-related solution set of which other solutions 
are special cases obtained by different choices of any 
permissible set of two of the remaining parameters. This 
solution is stated in theorem form as follows: 

Theorem 1 The Algebraic Riccati Equation 
for second-order mechanical systems is uniquely solved 
in closed form by the open-loop matrix pair 


where IT = — ©+ b^C+ a (2>+ Q*) 


Using the closed-loop formulation of (35), 
Theorem 1 can be alternatively stated as follows: 


( 36 ) 


Theorem 2 The Algebraic Riccati Equation 
for second-order mechanical systems is uniquely solved 
in closed form by the closed-loop matrix pair 
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(33) 

a M 



Q £ 0, K > 0 


(37) 


where FI = -0+ b^C + a(£> + 0. For K to be positive 
definite , the parameters of (33) are subject to the 
following constraints: 

1 

K B >0=>b>0; (K n -K B *KX)>0=>9>— * 

b (34) 

Note that the constraints (34) require a damping increase 
but, in general, allow positive, negative or zero stiffness 
change. 

The solution of Theorem 1 contains all possible 
sets of solutions to the optimal full-state-feedback control 
problem in which collocated perfect sensor-actuator pairs 
are present. This includes all combinations of linear 
structures and linear viscous dampers. The precise 
member (actuator-sensor pair) locations are set by the 
connection matrix The solution is best interpreted 
physically by regrouping terms in the closed-loop form of 
the original objective function as follows: 


1 


where ,0, (0 = K n )is an arbitrary matrix and a and b 
are arbitrary scalar constants. For K to be positive 
definite, the parameters of (35) are subject to the 
constraints (34), 


For stability, (9) requires Q, ^ 0 which translates, using 
the approach of (34), to 

a > 0 and (D >—9T 
b 

a very conservative requirement unless a is small relative 
to b. The above results are collectively summarized in the 
following theorem: 


V = — J’(z * Qz + u * Ru)dt = 
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Theorem 3 All control laws which can be 
implemented by collocated actuator-sensor pairs, 
including structural stiffness or linear viscous dampers, 
to minimize 

V = — J~ (z * Qz + u*Ru)dt, Q £ 0, R>0 
2 

subject to the equations of motion 


kit 


V 

UJ 


= -C(q q>‘ 


2a3C 2n 

2T\* 2bi>-2-dM 


Mi 


(35) 


^q+(2) + ^)q+^Cq = fti, z = 

are given by 

u = -aR '® * q + bR ’-S * q 

resulting in 

• • 

*Bu = -a#R * q - b'ER 1 # * q = A3(q + M>q 

where a and b are scalar constants . 

The closed-loop objective function is 
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v = 



A 



n 


2b©- 2ai( 



(38) 


A AAA 

where FI = - 0+ b3C+ a(©+ Q) and the value ofV is 


2aX a g 

ag 2b© - laid _ 

a© + bsc aid 
aid b id 


(40) 


V = — z * (0)Kz(0) 

2 

fq* (0)1 

U*(°)J (3 9) 

where 0 is a symmetric , positive-definite matrix . The 

A 

parameters a, b, and 0 are arbitrary but subject to the 
following constraints for a valid (i.e., K>0) and stable 

A 

(i.e., Q, > 0 ^solution: 


-(q - (0) q* (0))| 


0 a 9d 
a id b if 


* a * * a A 

b > 0 and 0 > — Od ; a > 0 and *D > ; 

b b 

The integrand of (38) is a "generalized-energy” 
function for linear, viscously-damped, multi-variable 
mechanical systems which are altered by changing member 

A A A 

damping and stiffness. For 0 = b3C , a=l/2, and © = 0, 
(3) reduces to the standard Lagrangian for conservative 
systems. However, since b>0 is required for a solution to 
exist, the closed-loop damping cannot be zero and a direct 
extrapolation to the Lagrangian is not valid. The 
expression for the control forces or system changes, 

• • 

(Su = -a«R * q - b«R~'s * q = A3Cq + A© q) indicates 
clearly that, in the diagonal gain format of Section I, the 
diagonal values of R are the inverse of the spring and 
damping constant values. Hence, a heavy penalty on 
control, R, results in a small control effect 

Physical Interpretation of General Solution 

Theorem 2 provides significant insight into the 
physical interpretation of possible optimal control 
solutions for structural systems. It gives key to the 
closed-loop (redesigned system) dynamics relationships 
between what is being minimized and the final value of the 
objective function in terms of mass, stiffness and damping 
of the final system. Of course, the final control law (i.e., 
the structure design change) still depends on the choice of 
R, < B, a, and b. By imposing additional constraints to the 

A 

equation set, other recognizable forms of Qand K can be 
obtained. Examples are as follows: 

AAA 

Example 1: Imposed constraint- 0 = a©+b3C. 
Solution: 


This solution is equivalent to the that of Koiitani and 
Ikeda. presented in open-loop form in (24) and (25). 

A A 

Example 2: Imposed constraints- 0 = bK, and a =0. 
Solution: 


x> 

II 
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o 

o 

1 

*> 

► o 
1 

[0 2b ©J 

.0 


(41) 


In this case, the integral of a constant times the Ravleieh 
Dissipation Function is minimized and the objective 
function cost is b times the total initial energy. 


A A 

Example 3: Imposed constraints-0 = bK 


Q = 2a 


K = b 


X -(©+ Q) 

2 

1a b A 

-(!>+£)* 

12 ^ a 

* a A 
X ~9d 
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a- 


-M 9d 
Lb 

For positive definite K, (38), the requirements are 


b > 0 and 


— I < 9d* X 
b 


(42) 

(43) 


Orthogonal modal vectors, <j» , may be chosen such that 
the mass and stiffness matrices are normalized to yield 

A A A A A A ^ a 

<t>* M d* = I and OflCd) - (Q) whcre q j s a diagonal 


matrix with diagonal elements co t the undamped natural 

frequencies of the closed-loop system. Substituting these 
orthogonality properties into (39) yields the condition for 
positive definite K to be 


b > 0 and X > 



(44) 


where ®* is the smallest eigenvalue of the undamped 
system. 


Consider the single-degree-of-freedom subset of 
this solution. In order for K to be positive definite, 
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|a| t 

— <»• = 11 T =H K , 


where the lower case symbols represent the single degree- 
of-freedom version of their upper case counterparts. The 
value of Kj 2 at which this occurs is independent of both a 
and b. K12 is constrained to be less than twice critical 
damping and (ial/b) is constrained to be smaller than the 
natural frequency. For a=0, the cost is the initial energy in 
the open-loop system. (The closed-loop system is damped 
even though damping does not appear in the K matrix.) 

A A 

Example 4: Imposed constraints- 9 = bK and 

A A 

© = (2a / b) 04 ■ Solution: 


Example 7: The solution of Belvin and Park . (29) 
through (32), is not a member of this solution set in 
general because it assumes values for all four terms of Qjj 
and for r. The resulting K22 is not necessarily 
constrained to the same predefined connection form 
although K12 is. Equivalence occurs in the following 
special case: 

(2a + b) , - „ 

* = - — - — -04, tj = a, 9 = 2b5C 

b (50) 

Note that the original statement of the Belvin -Park 
solution required to exist and T = 


Q - 2a 


-(*>+ Q) 
2 


-(*>+ Q*) 
.2 


In this case, if there are no gyroscopic terms, 

A A 

K = (b / 2a)Q and from equation (8), 

V = - z * (0)Kz(0) = — z* (0)Qz(0) = - J 0 ’z * Qzdt 
2 4a 2 

(47) 

This simple scalar relationship, which requires a>0, is 
interesting in that it is equivalent to the integrand 
remaining constant for a period of time equal to b/2a and 
then dropping to zero, or dropping linearly to zero over 
time period b/a. 

AAA 

Example 5: Imposed constraints- 0 = a2?+ b^C 
. * ( 2aY. 


04 -04 

b 


and © = ^ — Jo4. Solution: 


Q = 2a* ? K=b 


0 04 


In this case, the integral of the total energy is minimized . 

A 

Example 6: Imposed constraints- 0 arbitrary and a=0 


(-0+b^) 


-0+ba: 


A 

1 

„ L 0 

0 

K = b 


Lo 

b04_ 


In this case, note that if b = 1 and 0=1, the final cost, 
which is minimal for any particular $ and R, is the sum 
of the kinetic energy and the squared displacements. This 
has the potential of being a highly useful physical case 
where practical solutions exist. 


In the general LQR control design problem, 
choosing the control weighting matrix R (or r if the 
connection matrix $ is not predefined) is a difficult task. 
In the closed-form solution of Theorem 3, a permissible 

A 

set of (Q > 0) and/or 0 matrices can be selected based 

on the physics of the problem, and the equations solved for 
T, the optimal feedback gain to minimize the resulting 
function. Solutions, where they exist, indicate what the 
ideal control would be but are not always physically 
practical. The following subsections illustrate this 
process. 


Rewriting the Q equations of (37) gives insight 
into selection options for final solutions in terms of 
frequently used weighting matrix forms, where they exist 

Qn=2aflC + a J r (51) 

A 

Q l2 = - 0+ bX + a (2> + Q) + abr ^52) 

Q u = 2b *D + b 3 r — 2a 

Having already chosen two of the four parameter matrices 
allowable for a unique solution, only two remain. From 
(51) to (53) it is clear that one choice is to choose either 

A 

Qll or Q22* which fixes r, and either Q12 or 0(i.e., 
K11), either of which fixes the other. (Note that Q12 
must contain a term which equals the gyroscopic term 

a Q, if it exists, since all other matrices are symmetric and 
T is defined to be symmetric.) If the available 
connectivity, #, is known, then the optimal member 
gain is R‘* and is solvable from F. A second possibility 

A 

it to choose both Q 12 and © which fixes r and, hence, 
Qll and Q22- The difficulty in doing this is that the 
closed-loop objective function is of the generalized energy 
form and, consequently, many frequently used non-energy 
Q forms produce impractical solutions. 


Choosing both Qi i and Q22 assumes that an 
additional constraint equation, obtained by eliminating r 
between (51) and (53), is satisfied. This constraint, 

Q„=^(Q»-2bD + 2a*0 + 2a* 

b (54) 

implies a fixed relationship between the mass, stiffness, 
and damping matrices which is not likely to exist initially 
and may not be readily imposed on the system. 

Some examples, based on various choices of 
parameter matrices are presented below: 

Example 8: Qi 1 = I, Q12 = 0 
a 

b 2 , . 

Q m = 2b© - 2a^f + —(I - 2a»C) 
a 

* b T 

0 = -b9C + a© + — I 
a 

This solution is easily shown to be physically impractical 
by examining the closed-loop stiffness matrix. 

»C + A3C = -I-3C 
a 

For most problems, this matrix will require members with 
negative stiffness. For example, for the three-mass 
system, Configuration 1, in Fig. 3, referring to the 
notation of Fig. 2 and assuming : 

G k =[2,2,4, 1,3, 1] and a=0.2, 

the closed-loop spring gains become 

G k =[3, 3, 1,-1, -3,-1] 

with negative spring constants on the ungrounded springs. 
For the choice of Qn=I, R cannot be positive definite 
unless the closed=loop stiffness matrix is diagonal, that is, 
no masses are interconnected. 

A 

Example 9: Q u =a£, 0 = b9C 

1 

r = — © 

b 

a 2 

Q„ = © + 2aflC 

b 

Q h = b© - IzM 

In this case, r, (i.e., R) is negative and the stiffness and 
damping changes arc 


a 

A3C = — ©, and A© = -© 
b 

The closed-loop damping is zero and the system is not 
asymptotically stable. 

Example 10 : z*Qz = y*y = z*C*Cz 
(C=[ Cj c 2 ]*), arbitrary Q. 

This choice fixes all three Q matrices and, hence, imposes 
a requirement that constraint (54) be satisfied, that is 

a 2 a 2 

C, * C, + — C. * C = —(-2b© + 2a M ) + 2a9C 
' b 2 b 2 

Furthermore, (52) becomes 

r = — [C, * C, + 0- boa - a(© + $)] 
ab 

which requires Cj*C 2 to contain a skew symmetric term 

to offset the gyroscopic term. Even if the gyroscopic term 
is zero, the constraint between Cj and C 2 will be difficult 

if not impossible to meet in order to get a valid solution. 

Optimal R for Spe cified Closed-Loon Q 

Prespecifying closed-loop parameters in order to 
solve for the optimal R is an alternative to the open-loop 
case above. From (35), the appropriate relationships are 

Q„ = 2aac = 2aac + 2a 2 r (55) 

Q n =-e+a(©+$) + b«: 

= -0+a(© + £)+b3C + 2abr (56) 

Q n = 2b©— 2a!*f = 2b© — 2a^f+ 2b J T (5 7) 

Example 11 : Q„ = /, Arbitrary second choice. 

In (his case, we are minimizing the square of the 
closed-loop displacements. The closed-loop system 
requires 


That is, the closed-loop stiffness is diagonal, no matter 
what its open-loop form. Clearly, from the diagonal gain 
conceptualization of Section I, the solution is decentralized 
control with all masses independently grounded (and all 
interconnections eliminated). However, this is rarely 
practical in passive systems, if ever. In active systems, 
considerable energy must be expended to cancel the internal 
forces of the connecting springs and dampers, an inefficient 
approach at best. Of course, removing the 


interconnections and applying independent, decentralized, 
active controls (i.e. "flying in formation") will work but is 
also likely to be impractical. 


YL Concluding Remarks 

The design of dissipative structure/control 
systems using linear quadratic regulator (LQR) theory is 
investigated. Results are summarized as follows: 

1) The problem is formulated in physical 
coordinates, in term of mass, stiffness, and damping 
matrices. A diagonal-gain structural connectivity form is 
used to clarify physical interpretation. 

2) Solutions to the algebraic Riccati equation are 
developed which constrain system changes to predefined,, 
physically possible connections in mechanical systems. 

3) A general closed-form solution is presented 
which can be adapted to many LQR objective functions. 
The closed-loop objective function is a particular 
generalized energy form, regardless of what open-loop 
weighting matrix (Q) form is chosen. 

4) An optimal control gain weighting matrix (R) 
is shown to be available for properly specified objective 
functions. Improper choices of Q weighting matrices are 
shown to result in impractical hidden constraints or 
physically non-existent solutions for control gains. 


7. Belvin, W. K.; and Park, K. C.: Structural Tailoring 
and Feedback Control Synthesis: An Interdisciplinary 
Approach. AIAA Journal of Guidance, Control, and 
Dynamics, Vol. 13, No. 3, pp. 424-May-June, 1990. 

8. Inman, D. J.: Critical Damping in Complex 
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